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1. Introduction 

This paper is devoted to the simultaneous identification of one coefficient and the initial 
conditions in a reaction-diffusion system using the least number of observations as 
possible. 

Let C 1" be a bounded domain of M n with n < 3, (the assumption n < 3 is necessary 
in order to obtain the appropriate regularity for the solution using classical Sobolev 
embedding, see |3J. We denote by v the outward unit normal to Q on T = dQ assumed 
to be of class C l . Let T > and to G (0,T). We shall use the following notations 
Q = Q x (0,T), Q = Q x (t ,T), S = T x (t ,T) and S = T x (0,T). We consider the 
following reaction-diffusion system: 

d t u = Au + a{x)u + b(x)v in Q , 

d t v = Av + c(x)u + d(x)v in Q , . . 

u(t, x) = g(t, x), v(t, x) = h(t, x) on E , 

u(0,x)=uo and v(0, x) — Vq in Q. 

Reaction-diffusion systems are frequently used to model several physical applications, for 
example : in biology and medecine, emergence and growth of cancer and angiogenesis 
(see p[]), in ecology, prey-predator systems, insect dispersal (see [T%]). in chemistry, 
reaction in the presence of diffusion could produce spatial pattern of the chemical 
concentration (see [21J). 
Our problem can be stated as follows: 

Is it possible to determine the coefficient b(x) and the initial conditions uq(x), vo(x) for 
i6ll from the following measurements: 

t + T 

d t v\(t ,T)xuj and Au(T', •), u(T', •), v(T', •) in Q for T" = — - — , 
where w be a sub domain of Q ? 

Throughout this paper, let us consider the following set 

A(R) = {$ G L°°(fi); ||$|Uao (n) ^ R}, 
where R is a given positive constant. 

If we assume that (uo, Vq) belongs to (H 2 (Q)) 2 and g, h are sufficiently regular (e.g. 
3 e > such that g,h G H\to, T, H 2+E (dQ)) n H 2 (t , T, H £ (dtt))), then © admits a 
solution in H 1 ^^, H 2 (Q)) (see |Ej). We will later use this regularity result. 
Let (u,v) (resp. (m, u)) be solution of (JTJ) associated to (a, 6, c, cf, u , v ) (resp. (a, 6, 
c, d, Mo 5 ^o)) satisfying some regularity and "positivity" properties: 

a, b, c, d, 6 G A(i2), 

There exist two constants r > 0, c > such that 
< ~ ~ ( 2 ) 

u o > 0, m > r, c > c , 6 > 0, c + dr > 0, 

(7 > and h > r. 

Such assumptions allows us to state that the function v satisfies \v(x, T')\ > r > in fl 
(see H thm 14.7 p. 200). 
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We assume that we can measure d t v on uo in the time interval (t , T) for some to £ (0, T) 
and Am, u and v in O at time T" G (t ,T). Our main results are 

• A stability result for the coefficient 6(x) (or a(x)): 
For uq, vo in H 2 (Q) there exists a constant 

(7 = C{Q,u,c ,t ,T,r,R) > 

such that 

\b - b\ 2 L2{n) < C\d t v - d t v\ 2 L2{[t0>T)xij) + C\Au(T', •) - Au(T', .)\ 2 L2(n) 
+C\u(T', •) - u(T', .)\l Hu) + C\v(T', ■) - v(T', -)|i 2(n) . 

• A stability estimate for the initial conditions uq, vq\ 
For wo, vo, uq, vo in H 4 (Q) there exists a constant 

C = C{Q,u,c ,t ,T,r,R) > 

such that 

where E = - d t v\h {{ t ,T)xu>) + -)|^ (n) + |^(T', ■) — v(T, .)|^ (n) . 

The key ingredient to these stability results is a global Carleman estimate for a two 
by two system with one observation. Controllability for such parabolic systems has 
been studied in pQ. The Carleman estimate obtained in [I] cannot be used to solve 
the inverse problem of identification of one coefficient and initial conditions because 
of the weight functions which are different in the left and right hand side of their 
estimate. We establish a new Carleman estimate with one observation involving the 
same weight function in the left and right hand side. We prove a stability result for the 
initial conditions following the method of [22] • Concerning the stability of the initial 
conditions we use an extension of the logarithmic convexity method (see ^D])- We can 
also cite he provides an estimate for the initial condition for a general parabolic 
operator, while the logarithmic convexity method works only for self-adjoint operators 
(see section 2 and Theorem 3 in |15j). 

The simultaneous reconstruction of one coefficient and initial conditions from the 
measurement of one solution v over (to,T) x uj and some measurement at fixed time T' 
is an essential aspect of our result. In the perspective of numerical reconstruction, such 
problems are ill-posed. Stability results are thus of importance. 

For the first time, the method of Carleman estimates was introduced in the field of 
inverse problems in the work of Bukhgeim and Klibanov [Sj; also see, e.g., jlj, ^3] and 
[Tl] for some follow up publications of these authors. So far, the method of j3] is the 
only one enabling to prove uniqueness and stability results for inverse problems with 
single measurement data in the n-dimensional case with n > 2, which has generated 
many publications, including this one. While this method can provide Holder stability 
results, the topic of the Lipschitz stability is a more delicate one. The first Lipschitz 
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stability result for a multidimensional inverse problem (for a hyperbolic equation) was 
obtained by Puel and Yamamoto [19J, using a modification of the idea of |5,. In JJJ 
this result was extended for a parabolic equation. The main difference between our 
work and [UJ is that we consider a coupled system of parabolic equations, and the 
additional data are given only for one component of this system, the function dtv(x,t) 
for (x,t) G lj x (t , T) . Inverse problems for parabolic equations are well studied (see 
El, [TT] . |22j). A recent book of Klibanov and Timonov 16J is devoted to the Carleman 
estimates applied to inverse coefficient problems. In our knowledge, there is no work 
about inverse problems for coupled parabolic systems. 

The used method allows us to give a stability result for the coefficient a{x) adapting 
assumption 13.11 On the other hand, since we only measure d t v on lj, we cannot obtain 
such stability results for the coefficients c(x) or d(x) of the second equation of (JTJ). 
For the reconstruction of two coefficients the problem is more complicated. We obtain 
partial results with restrictive assumptions on the coefficients a(x), b(x), c(x) and d(x). 
In order to avoid such assumptions, we think it is necessary to use other methods such 
as those used in [T2j . 

Our paper is organized as follows. In Section 2, we derive a global Carleman estimate for 
system (JTJ with one observation, i.e. the measurement of one solution v over (t , T) x lj. 
In Section 3, we prove a stability result for the coefficient b(x) when one of the solutions 
v is in a particular class of solutions with some regularity and "positivity" properties. 
In Section 4, we prove a stability result for the initial conditions. 

2. Carleman estimate 

We prove here a Carleman-type estimate with a single observation acting on a subdomain 
lj of Q in the right-hand side of the estimate. Let us introduce the following notations: 
let lj' (£ uj and let (3 be a C 2 (Q) function such that 



P > 0, in Q, (3 = on dtt, mm{\V(3(x)\, i£fl \ lj'} > and d u (3 < on dtt. 

Then, we define /3 = (3 + K with K = m||/3||oo and m > 1. For A > and t G {t ,T), 
we define the following weight functions 



<p(x,t) 



e Xf3{x) 



T](x,t) 



e 2\K _ e Xf3(x) 



(t-t )(T-t) 



(t- to )(T-ty 



If we set ip = e sv q, we also introduce the following operators 
MiV>= - Aifj - s 2 \ 2 \V(3\ 2 (p 2 ilj + s(d t r])ip } 
M 2 ip = d t ip + 2s\ipV(3.V^ + 2sAV|V/3|V- 

Then the following result holds (see jS]). 
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dt 



(3) 



Theorem 2.1 There exist X = X (Q,uj) > 1, Sq — s o{^o,T) > 1 and a positive 
constant Cq = Co(Q,u,T) such that, for any A > Ao and any s > s , the following 
inequality holds: 

11^-^)11^ + ||M 2 (e-^)||l 2(Q) 
-sA 2 J J e- 2s \\Vq\ 2 dx dt + s 3 \ 4 J J e~ 2s Vk| 2 dx dt 

s 3 A 4 j j e" 2s Vk| 2 dxdt+ jj e~ 2sri \d t q - Aq\ 2 dx 

J to J u) J J Q 

for all q G H\t , T, H 2 (ty) with q = on S. 

From the above theorem we have also the following result (see jU] and [TT] ) 



Proposition 2.2 There exist X = \ (Q,u>) > 1, Sq — Sq(Xq,T) > 1 and a positive 
constant Cq = Co(Q,oo,T) such that, for any A > Ao and any s > sq, the following 
inequality holds: 



e- 2s \-\\d t q\ 2 + \Aq\ 2 ) dx dt 



+ 



sA 2 f j e- 2s \\Vq\ 2 dx dt + s 3 X 4 j j e~ 2s \ 3 \q\ 2 dx dt 

JJq JJq 

s 3 X 4 ft e- 2s \ 3 \q\ 2 dxdt+ jj e~ 2sr) \d t q - Aq\ 2 dx 
JtoJuj JJq 



<c 

for all q G H l {t , T, H 2 (J})) with q = on E. 

We consider the solutions (u, v) and (u, v) to the following systems 

d t u = Au + au + bv in Q , 

d t v = Av + cu + dv in Q , 

u(t,x) — g(t,x), v(t,x) — h(t,x) on S , 

u(0,x) = u and v(0,x)=v in Q, 

d t u = Au + au + bv in Q , 

d t v = Av + cu + dv in Q , 



dt 



and 



u(t,x) = g(t,x), v(t,x) = h(t,x) on S , 
u(0,x) = u and v(0,x) = v in Q. 



(4) 



(5) 



(6) 



We set U = u — u, V = v — v , y = d t (u — u) , z = d t (v — v) and 7 = 6 — 6. Then (y, z) 
is solution to the following problem 

d t y = Ay + ay + bz + jd t v in Q , 

d t z = Az + cy + dz in Q , 

y(t, x) — z(t, x) — on S , (7) 

?/(0, a;) = AU(0,x) + aU{0,x) + bV(0,x) +jv(0,x), in Q, 

^ z(0,x) = AV(0,x) +cU{0,x) +dV(0,x) in ft. 
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Note that the previous initial conditions are available for all T' G (0, T). Indeed with 
(JHJ) and © we can determine y(T',x), z(T',x) and we obtain 

y(T', x) = AU(T', x) + aU(T', x) + bV(T', x) + ^v(T', x), 

z(T', x) = AV(T', x) + cU(T', x) + dV(T', x). 
We consider the functional 

I(q) = s- 1 II e~ 2s ^-\\d t q\ 2 + |Ag| 2 ) dx dt 



+s\ 2 / / e~ 2s W Vg| 2 dx dt + s 3 A 4 / / e~ 2sr, ^\q\ 2 dx dt. 
JJq JJq 

Then using the Carleman estimate (J3J), the solution (y, z) of (|IJ) satisfies 

i{y) + i(z) < Ci[s 3 A 4 I I e - 2s VM 2 dx dt 

JtnJuJ 



'to Jul 

+s 3 A 4 / / e" 2s VM 2 dx dt 

JtoJuJ 



+ JJ e- 2sr > (\ay\ 2 + \bz\ 2 + \-fd t v\ 2 ) dx dt + JJ e" 2 ^ (|q/| 2 + \dz\ 2 ) dx dt] 
Let £ be a smooth cut-off function satisfying 



£(x) = 1 ViG u/, 
< £0) < 1 Vx G to", 
£0) = Vx 6 R" \ a/', 



where w' g w" g w g O. 

We shall estimate the following three terms 



I := s 3 A 4 I I e~ 2s \ 3 \y\ 2 dx dt, 

J := II e~ 2sri \bz\ 2 dx dt or /Y e" 2 ^ \dz\ 2 dx dt, 
JJq JJq 

K := II e~ 2sv \ay\ 2 dx dt or II e~ 2sr] \cy\ 2 dx dt. 



JJQ JJQ 

The main difficulty is in estimating from the above the term /, because values of the 
function y(x,t) are unknown for (x,t) G to x (to,T). In doing so, we rely on the 
assumption c(x) > Co where Co is a positive constant, see (jHJ). For the first term I, we 
multiply the second equation of (J7J) by s 3 A 4 e~ 2sT? £<^ 3 2/ and we integrate over (t ,T) x to. 
We obtain 

I' := s 3 A 4 I I ce- 2sv ^ 3 \y\ 2 dx dt = s 3 X 4 I I e" 2 ^ 3 ^ - Az - dz)y dx dt 

J ZqJ UJ J IqJ uj 
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s 3 A 4 / / e- 2sr >& 3 (d t z)y dx dt - s 3 A 4 / / e~' 2sr ' '(Az)y dx dt 



to Jul 




t Jul 



-s 3 A 4 / / de-'^&zy dx dt = h + I 2 + h 




to Jul 



By integration by parts with respect to the time variable, the first integral, ii, can be 
written as 



h = s 3 X 4 



j j e- 2si >& 3 z(d t y) dx dt + 2s 4 X 4 f f e~ 2sT1 ^ 3 {d t r]) zy dx dt 

J to Jul J to Jul 



3s 3 A 4 / / e- 2sr >& 2 (d t (p)zy dx dt. 

J to Jul 



We write Ii = I{ + If with 



I\ = -s 3 X 4 / / e~ 2sri & 3 z{d t y) dx dt, 

JtoJul 

II = 2s 4 A 4 I j e- 2sri ^ 3 (d t ri) zy dx dt - 3s 3 A 4 / I e- 2sr >^ 2 (d t ip) zy dx dt. 

J t()J UJ J t()J UJ 

Using Young inequality, we estimate the two integrals I\ and If. We have 



\ll\ <s 3 \ 4 



< -s 7 X 8 
- 2e 



C £ s 4 X 4 ft e- 2s ^V\z\ 2 dx dt + es~ 4 \- 4 fl e~ 2s ^ ^ x \d t y\ 2 dx dt 

JtoJul JtoJul 




e 2s VH 2 dx dt + 



to Jul 



e- 2sr '^- l \d t y\ 2 dx dt. 



Q 



The last term of the previous inequality can be "absorbed" by the terms in I(y) for e 
sufficiently small (for exemple e — |). 

fT 



4x4 



|/f I < Cs 4 A 




+s~ L X 



i \-i 




sX I / e- 2s ^ 2 ip(if 2 \d t r]\ 2 + \d t <p\ 2 )\z\ 2 dx dt 

toJ ui 



e- 2sTI V 3 \y\ 2 dx dt 



to Jul 



< C 



s 5 X 5 / / e" 2s VM 2 dx dt + s 3 A 3 / / e~ 2s \ 3 \y\ 2 dx dt 



t J Ul 



The last inequality holds through the following estimates 

\dt<p\ < c(n, uj)Tip 2 , \d t r]\ < c(n, uj)t<p 2 , if < c(n, u)t 4 ^ 3 . 

The last term of the previous inequality can be "absorbed" by the terms in I(y) for s 
and A sufficiently large. Finally, we obtain 

< Cs 7 A 8 I j e- 2sv cp 7 \z\ 2 dx dt + "absorbed terms" , 

J to J Ul 



Inverse problems for a reaction- diffusion system 



8 



where C is a generic constant which depends on Q, u and T. 

Integrating by parts the second integral I 2 with respect to the space variable, we obtain 



J 2 = - S 3 A 4 f j A(e- 2si >£ip 3 y)z dx dt. 

J to J OJ 

If we denote by P = e _2s,? £<^ 3 , then we have 

T 

I 2 = - S 3 A 4 j j (PAy + 2VPVy + yAP)z dx dt. 
We compute VP and AP and we obtain the following estimation for I 2 



h\ < 5 3 A 4 



es- 4 \^ 

T 




e-^cp-'lAy] 2 dx dt + C £ s 4 \ 4 / / e~ 2s \ 7 \z\ 2 dx 



U)J Ul 




+es- 2 \- 2 e~ 2s \\Vy\ 2 dx dt + C £ s 2 \ 2 e - 2s \ 5 \z\ 2 dx dt 



t(,J Ul 




IT 



+e / / e- 2s VM 2 dxdt + C £ / e~ 2s Vbl 2 dx dt 



Therefore we obtain 



II 



e-^ip' 1 1 Ay | 2 dx dt + s\ 2 / / e~ 2sri V \Vy\ 2 dx dt 



Q 



Q 



->-\V j j e~ 2s VM 2 dx dt 



a 



s 7 X 8 I I e-*VM" dx dt 




-2sr) ir 7\\2 



to Jul 



IT 

J toJ ul 



+s 5 A 6 / / e- 2s \ 5 \z\ 2 dx dt + s 3 X 4 / / e" 2s VM 2 dx dt 




t Jul 



The first three integrals of the r.h.s. of the previous inequality can be "absorbed" by 
the terms in I(y) for e sufficiently small. Finally, we have 

\h\< Cs 7 \ 8 j j e~ 2s \ 7 \z\ 2 dx dt+ "absorbed terms" . 

J toJ Ul 

For the last integral I3, we have 



3\4 



|/ 3 1 < Cs 6 \ 



C £ e~ 2s \ 3 \z\ 2 dxdt + e // e" 2s VM 2 dx dt 



to Jul 



Q 



Finally, if we assume that there exists a constant Co > such that c > Co in u, we 
have thus obtained for A and s sufficiently large and e sufficiently small the following 
estimate: 



|/| < < £/a s 
Co c 




e" 2s VM 2 dx dt. 



(9) 



to J Ul 



For the integrals J and K, since a, b, c, d € A(_R) and using the estimate 

1 < C(Q, cu)TV/4, 
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we have 

\J\<C [[ e~ 2sv if 3 \z\ 2 dx dt, 



K\<C / / e- 2sr > <p 3 \y\ 2 dx dt, 



and these terms can be "absorbed" by the terms I(y) and I(z) for A and s sufficiently 
large. If we now come back to inequality (JEJ), using the estimates for I, J and K, and 
choosing A and s sufficiently large and e sufficiently small, we can thus write 

I{y) + I{z) < CMs 3 A 4 CI e- 2s \ 3 \z\ 2 dx dt + s 7 \ 8 f I e~ 2s VM 2 dx dt 

+ J J e- 2sr > | jd t v\ 2 dx dt}. 
Observing that 

if 3 < c(n,oo)T\ 7 , 

We have thus obtained the fundamental result 

Theorem 2.3 We assume a, b, c, d G A(R) and that exists Cq > such that c > Cq 
in to. Then there exist \\ = Xi(Q,u) > 1, s\ — si(X\,T) > 1 and a positive constant 
C\ = Cx(Q, D, c , R, T) such that, for any A > Ai and any s > s\, the following inequality 
holds: 

I{y)+I{z) < d s 7 \ s [ [ e" 2s VM 2 dx dt+ {I e- 2sTI \-fd t v\ 2 dx dt ,(10) 

Jt Juj JJq 

for any solution (y, z) of (0). 

3. Uniqueness and stability estimate with one observation 

In this section, we establish, a stability inequality and deduce a uniqueness result for 
the coefficient b. This inequality (|T5|) estimates the difference between the coefficients 
b and b with an upper bound given by some Sobolev norms of the difference between 
the solutions v, and v of (J3j) and (JBJ). Recall that U = u — u, V = v — v, y = d t {u — u), 
z = d t {v — v), 7 = 6 — 6 and 

d t y = Ay + ay + bz + ~fd t v in Q , 

d t z = Az + cy + dz in Q , 

y(t, x) = zit, x) = on S , 

y{0,x) = AU{0,x) +aU{0,x) +bV{0,x) +jv{0,x), in Q, 

k z(0,x) = AV(0,x) +cU(0,x) +dV(0,x) in a 

The Carleman estimate (jlO)) proved in the previous section will be the key ingredient in 
the proof of such a stability estimate. 

Let T' — |(T + to) the point for which $(t) = ^ t _ to ^ T _ t ^ has its minimum value. 
For {u, v) solutions of Q, we make the following assumption: 
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Assumption 3.1 There exist r > 0, Cq > such that b > 0, c > c , c + dr > 0, Uq > 
0, Wo > r, g > and h > r. 

Such assumption allows us to state that the solution v is such that \v(x,T')\ > r > 
in Q (see [20]; theorem 14.7 p. 200). Furthermore if we assume that Uq, Vq in H 2 (Q),the 
solutions of (jHJ) belong to -ff 1 (to, T, H 2 (Q)). Then using classical Sobolev embedding 
(see [3J), we can write for n < 3, that c^v belongs to L 2 (to, T, L°°(0)) and we assume 
that |<9^| i2(t0)T) g A(-R). 
We set ^ = e~ sv y. With the operator 

M 2 ^ = fltV + 2sA(^V/?.V^ + 2s\ 2 ip\V/3\ 2 ip, (11) 

we introduce, following [2|, 



J = 3? / / M 2 ^ V dardt 
We have the following estimates. 

Lemma 3.2 Lei A > Ai and s > s% and let a, b, c, d G A(R). We assume that 
assumption \3. 1\ is satisfied then there exists a constant C = C(Q,u,T) such that 



111 < Cs- Z ' 2 \- 



s 7 X 8 I I e~ 2s VM 2 dxdt+ I f e- 2sr >\<y\ 2 \d t v\ 2 dx dt 



Proof: 

Observe that 



X\ < s-V 2 \- 2 ( f T f |M 2 ^| 2 dx dt) (s 3 X 4 f T f 
\Jt Jn I \ J t Jo. 



V2 / T , \ 1/2 

s 3 A 4 / / e- 2sr >\y\ 2 dx dt 
t Jn 



thus using Young inequality and the estimate 1 < C'T ip , we obtain 

\1\ < Cs-^ 2 \- 2 (\M 2 i,\ 2 L2{Q) + S 3 A 4 H e- 2s \ 3 \y\ 2 dxdt 
V " JtJ n 

which yields the result from Carleman estimate ([10|h 



Lemma 3.3 Let A > Ax, s > s\ and let a, b, c, d G A(R). Furthermore, we assume 
that uq, Vq in H 2 {VL) and the assumption ^. 1\ is satisfied. Then there exists a constant 
C = C(Q,u,T) such that 



| 7 £(T',x)| 2 dx 



(12) 



< Cs- 3 / 2 \- 2 



S 7 X 8 I / e -2 S ^7| z |2 dx dt 

trJ LD 



[ [ e- 2sr, \-f\ 2 \d t v\ 2 dx dt 



-C [ e - 2sv( - T '' x) | A U (T', x) + aU(T', x) + bV(T', x)\ 2 dx. 
Jq 



Inverse problems for a reaction- diffusion system 11 
Proof: 

We evaluate integral 1 using (fTTj) 

T= \ I I d t \i)\ 2 dxdt + s\ [ [ ipVP- V\ip\ 2 dxdt + 2sX 2 [ [ f3\ 2 \i)\ 2 dxdt 



2 J^J-i- „\ / / V7 f ,^\~7 Q\\ n i.\2 J™J4- i Q~\2 / / ,„\T7 Q\2\ n L\2, 



2 

= 1 I I W dxdt - sX I I V ■ {cpVP)\ip\ 2 dxdt + 2sX 2 I I (p\VP\ 2 \ip\' 2 dxdt, 
2 Jt Jn Jt Jn Jt Jn 

by integration by parts. With an integration by parts w.r.t. t in the first integral, we 

then obtain 

\ [ \^(T',.)\ 2 dx = l-sX 2 [ '[ if\V(3\ 2 \ij\ 2 dxdt + sX [ '[ (p(A/3)\ip\ 2 dxdt 
2 Jn Jt Jn Jt Jn 

since ip(t ) = and V<£> = 

Then, we have 

I ( e -^)(T>)| y ( T / ):c )|2 dx < 2 |J|+CsA(A+l) [ T [ e- 2s ^' x) ip\y\ 2 'dxdt. (13) 

Using Lp < ^-tp 3 the last term in (J13)) is overestimated by the left hand side of (|l(Jj) 
and this last one is absorbed by the l.h.s. of the inequality obtained in lemma I3~21 
If we now observe that 

y(T', x) = AU(T', x) + aU(T', x) + bV(T', x) + yv{T' x), 

we have 

HT',x)\ 2 > ~\ 1 v(T',x)\ 2 -\AU(T',x)+aU(T',x) + bV(T',x)\ 2 . 

m 

The regularity of the solutions of (JHJ) allows us to write that for n < 3, d t v is an element 
of L 2 (t ,T,L°°(n)). So, from \v(x,T')\ > r > 0, we have 

3 k G L 2 (t ,T), \d t v(x,t)\ < k(t)\v(x,T')\, Vi 6 H, t e (t ,T). 

Hence f|12|) can be written 

(e-^XT'^M'I^TOI 2 cfe 

n 

T 



< Cs^X- 2 [ [ e-^M'IMt^feT')! 2 dxdt 
Jta Jn 



'to Jn 

+CV 1 / 2 A 6 H e~ 2s \ 7 \z\ 2 dx dt 



tiy CO 



+C [ (e- 2sri){r ' x \\AU(T',x)\ 2 + \U(T',x)\ 2 + \V(T' J x)\ 2 ) dx. 
Jn 



Since k 6 L 2 (t ,T) implies that / |&(t)| 2 dt < k < +oo. For A large enough, the 

Jt 

term (1 — Cs~ 3 / 2 A -2 A;o) can be made positive: 
1 - Cs- 3 / 2 X- 2 k > C 2 > 0. 
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Using the fact that e" 2 ^*^ < e" 2 ^'^ Vx efi, Vt G (* , T), we deduce that 

;i_ Cs - 3 / 2 A- 2 A;o) / e- 2s "( T '^| 7 | 2 ^ 
Jq 



r 2 ' 



< Cs ll / 2 X 6 / / e~ 2s Vkl 2 c/x dt 



+C [ e- 2sr >( T '> x \\AU(T',x)\ 2 + \U(T',x)\ 2 +\V(T',x)\ 2 )dx, 
Jq 

where we have also used that \v(x,T')\ > r > in Q. Then, by virtue of the properties 
satisfied by (p and 77, we finally obtain 

|7|^ ( ,)<^ U/2 A 6 jjjzf dxdt (14) 

+-£r [ (\AU(T',x)\ 2 + \U(T',x)\ 2 + \V(T',x)\ 2 )dx. 
r ^2 Jq 

With (JHJ), recalling that U = u — u, V = v — v, y = d t {u — u) and z = d t {v — v), we 
have thus obtained the following stability result. 

Theorem 3.4 Let u be a subdomain of an open set Q o/M n ; let a, b, c, d £ A(R). 
Furthermore, we assume that Uq, vq in H 2 (Q) and the assumption \S. 1\ is satisfied. Let 
(u, v), (u, v) be solutions to d3J)-(0J). Then there exists a constant C 

C = C{Q,u,c ,t ,T,r,R) > 

such that 

\b-b\ 2 L2{n) < C\d t v - d t v\ 2 LH(t0tT)xuj) + C\Au(T', ■) - Au(T' , .)& (n) (15) 
+C\u(T', •) - u(T', .)\ 2 L2{n) + C\v(T', •) - v(T', .)\ 2 L2{n) . 

Remark 3.5 If we assume that u{T',-) = u(T',-) and v(T',-) = v(T',-) (such an 
additional assumption is sometimes made, e.g. in }Hf)- then the stability estimate 
becomes 

\b-b\ 2 L2{n) < C\d t v - d t v\ 2 LH{t0jT)xuJ y 

With Theorem 13.41 we have the following uniqueness result 
Corollary 3.6 Under the same assumptions as in theorem and if 

(d t v — d t v)(t, x) = in (t , T) x u, 

Au(T', x) - Au(T', x) = in Vt, 
v(T', x) — v(T', x) — in Q, 

Then b = b. 
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In this section, we use the same method as in |22j to state a stability estimate for the 
initial conditions uq, Vq. The idea is to prove logarithmic-convexity inequality. The 
following method has been used to obtain continuous dependence inequalities in initial 
value problems. If (y, z) is solution of JZJ), we introduce (yi,zi) and (2/2,-22) that satisfy 

d t y x = Ay 1 + ay x + bz x + jd t v in Q , 
d t z\ = Azi + cyi + dz\ in Q , 



yi(t,x) = zi(t,x) = 
2/i(0, x) = 0, 
[ zi(0,z) = 



on S , 
in Q, 
in Q, 



(16) 



and 



in Q , 
in Q , 
on S , 



d t y 2 = Ay 2 + ay 2 + bz 2 
d t z 2 = Az 2 + cy 2 + dz 2 
y 2 (t,x) = z 2 (t,x) = 

2/ 2 (0, x) = AU(0, x) + aU(0, x) + bV{0, x) + 7^(0, x) in fi, 
k z 2 (0,x) = AV(0,x) + cU(0,x) \ + dV(0,x) in Q. 

Then, we have 

2/ = 2/i + 2/2 and z = z 1 + z 2 . 

In a first step, we give an L 2 estimate for (2/1, z\) 

Lemma 4.1 Let a, 6, c, <i, l^t^^T) £ A(i?). Then there exists a constant 



(17) 



C = C(t ,T',R) > 0, 



such that 



(19) 



l3/i(*)ll a( n) + \zi(t)\l, (a) < C| 7 |i 2(Q) , t < t < T'. 
Proof: 

We multiply the first (resp. the second) equation of (fT6j) by 2/1 (resp. by z\). Then, 
after integrations by parts with respect to the space variable, we obtain 

\d t [ + \ Zl \*) dx = - [ (|V2/ 1 | 2 + |Vz 1 | 2 )rfx 
1 Jn Jn 

Jn Jn Jn Jn 

We use Cauchy-Schwarz and Young inequalities and we integrate over (to,t) for to < 
t < T' and we obtain 



l?/i(^)li 2 (n) + l z i(^)li 2 (n) <C 2 |7l + Ci / (\yi(s)\h m 

t 



+ \zi(s)\ 2 L 2 {n) )ds. 



The result follows by a Gronwall inequality. 
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In a second step, we use a logarithmic-convexity inequality for (y2,z 2 ) 

Lemma 4.2 Let a, b, c, d G A(R) and u , v , u , vq in H A (Q). Then there exist 
constants M > 0, C = C{R) > and C 1 = Ci(t , T', R) > such that 

\y2(t)\h m + \z 2 (t)\l m < C.M^iMT')^ + \z 2 (T')\l 2{Q) )^\ (20) 

(e~ < ^*° e^^) 

for t < t < T', where jj(t) = - _ cto - _ c . 

\e e j ) 

Proof: 

The proof of this lemma is just an application of Theorem 3.1.3 in [HI]- In fact, system 
(fTTj) can be written in the following form 

d t W + AW = BW, in Q , 
W{t,x) = on S , 

W(0,x) = W (x) in ft, 



where 



The operator A is symetric and the solution W satisfies 

\\d t W + AW\\ L2{n) <a\\W\\ LHn) , 

where a = \ \B\\ L oo^ < +oo since a, b, c and d are in A(R). If we assume that u , v , 
Uo, Vq are in H 4 (Q), the hypothesis of Theorem 3.1.3 in 10, are satisfied, thus we have 

\\w(t)\\L^ } < CiII^^II^II^ctoii^, 

With M*) = (g-Cio-g-CT^ - 

Since G C(t ,T, L 2 (Q)), we have ||W(fo)IUa(n) < Afi, and the result follows. 

■ 

The two previous lemmas allow us to prove the following 

Theorem 4.3 Let u be a subdomain of an open set Q o/M n ; let a, b, c, d G A(R). 
Furthermore, we assume that Uq, Vq, uq, vq in H 4 (Q) and Assumption \S. 1\ is satisfied. 
Let (u, v), (u, v) be solutions to (Ej)-(G|). We set 

E = \d t v - d t v\ 2 L2{(tQjT)xuj) + \u(T', •) - u(T', -)\ 2 H2{n) + \v(T', •) - v(T', -)|^ (n) . 

Then there exists a constant C = C(VL, u, c , to, T, r, R) > such that 

C 

K ~ U \ 2 L 2 {n) + \V - V \ 2 L 2 {n) < hgE y f° r < E < 1 
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Proof: 

Since v G H\t ,T, H 2 {Q)), we have v e £,°°(fi). In view of (HHJ), inequalities fl3J, l|20|) 
imply 

|y(V)l <2(|j/i(t,-)| 



l7li»(n) + M^{\y 2 {T\ -)lia (n) + |* 2 (T', .)li» (n) )^ 



Now, with (fTHjl . we write yi — y — y\ and z 2 = 2 — 21. Inequality (fTUj) gives us an 
estimation of \yi(T', -)li 2 (r>) and \z\(T', -)lL 2 (n) m terms of M^^- Then the definition 
of y and z in gives us an estimation of \y(T', -)|^ 2 
\U(T, -)\ 2 H2m and |V(T', Ol^n)- Finally, we obtain 



and \z(T', Oliam) m terms of 



lv(V)l 



L 2 (n) 



7Il»(0) + ^(|7i!a (n) + |t/(T', .)|^ fn) + -)\ 2 H , 



(«) 



2(H), 



(a similar estimate is obtained for ^(t, -)| 2 ). Following |22j with a time translation, if 
we use (JI5|) . the last estimate yields 



l«0 - Wo|ia( fl ) + l«0 - ^o|| 2 (r>) = |E/(t , -)li 2 (Q) + l^(*0, OIL 



/ y(s,-)ds + U(T', 

Jto 



\L 2 {n) 



+ 



to 



< m 3 / ds + c 2 e < a 



to 



E- 1 
log£ 



+ C 4 E < 



z(s,-) ds + V(T',-)\l 2(n) 
C 



\ogE\ 
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